Enhancing the superconducting transition temperature of BaSi2 by structural tuning 
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We present a joint experimental and theoretical study of the superconducting phase of the lay- 
ered binary silicide BaSi2. Compared with the layered AIB2 structure of graphite or diboride-like 
superconductors, in the hexagonal structure of binary silicides the sp 3 arrangement of silicon atoms 
leads to corrugated sheets. Through a high-pressure synthesis procedure we are able to modify the 
buckling of these sheets, obtaining the enhancement of the superconducting transition temperature 
from 4K to 8.7 K when the silicon planes flatten out. By performing ab initio calculations based on 
density functional theory we explain how the electronic and phononic properties of the system are 
strongly affected by changes in the buckling. This mechanism is likely present in other intercalated 
layered superconductors, opening the way to the tuning of superconductivity through the control of 
internal structural parameters. 



Nowadays, an important part of the activity research 
on superconductivity is focused on intercalated layered 
crystal structures. In these systems, where the relevant 
features for superconductivity are intrinsic to the lay- 
ers, the bond buckling of the atoms forming the layers is 
a structural parameter known to damage the supercon- 
ducting properties, regardless of the nature of the pairing 
mechanism The case of flat boron sheets in MgE$2 is 
well understood: superconductivity arises from a strong 
coupling between the sp 2 a-bonding intralayer electrons 
of boron and the in-plane bond stretching phonons @- 
4]. Besides, a buckling of the boron honeycomb struc- 
ture, as observed in ReB 2 was proved to decrease 
the superconducting transition temperature T c . Never- 
theless, recent studies on graphite intercalated supercon- 
ductors, namely (Yb,Ca)Cg 0, and the ternary silicide 
CaAISi @,Hj], characterized by large T c s, point to the im- 
portance in these systems of the electron-phonon (EP) 
coupling between the intcrlayer electrons and the out- 
of-planc vibrational modes of the atoms composing the 
layer. In that case, the buckling phonon modes corre- 
sponding to the antiphase motion along the c-axis of the 
atoms in the sheets can lead to an enhancement of the 
EP coupling. Experimentally, the effect of buckling can 
either be explored by means of high-pressure or directly 
by chemistry. An excellent testbed for such studies is 
the trigonal phase of binary silicides, like CaSi2 or BaSi2, 
which constitutes a family of layered intercalated super- 
conductors closely related to the graphitc/diboridc inter- 
calated systems. Unlike the so-called AlB2-like structure 
of graphite/diboride compounds, where the planes are 
flat, in the hexagonal structure of binary silicides the Si 
planes buckle. 

In this Letter we investigate superconductivity in 
BaSi2 and, in particular, we focus on its layered struc- 
ture (EuGe2-type structure, P-3ml space group) [9L[l"cij]. 



As it is sketched in the inset of Fig. [2j trigonal BaSi2 
is made of planes of Ba arranged in a triangular lattice, 
interspersed with buckled hexagonal planes of Si. This 
phase is metallic and was found to be superconducting 
with a critical temperature of 6.8 K [ll[. There are all 
the reasons to expect that BaSi2 is a standard s-wave su- 
perconductor, where the mechanism of superconductivity 
can be obtained by studying the EP coupling. However, 
recent theoretical studies [12| of trigonal BaSi2 yielded a 
value of T c almost an order of magnitude smaller than 
the early experimental finding of 6.8 K. These studies 
were based on state-of-the-art linear-response calcula- 
tions within density- functional theory (DFT), that have 
proved to describe very accurately the superconducting 
properties of many other similar materials. [l3j ]. Anal- 
ogous calculations for CaSi2 in its trigonal high pres- 
sure phase gave also a significant underestimation of 
T c [14J . Interestingly, the application of pressure on trig- 
onal CaSi2 was shown to increase the superconducting 
temperature up to 14 K [l5j. This finding was discussed 
in terms of a structural phase transitions under pressure, 
which is expected to soften the phonon modes associated 
to the collapse of the corrugated Si planes [l6| . 

In order to elucidate the origin of superconductivity in 
binary silicides, and at one time to shed light on the dis- 
crepant results found in literature, we performed a joint 
experimental and ab initio study of the superconduct- 
ing phase of BaSi2- Our experimental data show that 
T c varies strongly with the degree of buckling of the Si 



planes. As in the case of CaSi2 [15|, LUj], by flattening 
the Si planes, we increase the superconducting transi- 
tion temperature T c from 4 K up to 8.7 K. Our extensive 
DFT calculations allow us to determine how the buck- 
ling of the Si planes impacts the electronic and phonon 
band structures, and finally the EP coupling. Through 
the comparison of experimental and theoretical results 
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FIG. 1: (Color online) Left panel: Temperature dependence of 
the ZFC dc-magnetic susceptibility measured at 20 Oe of two 
trigonal samples of BaSi2 marked by their z coordinate (see 
text). For convenience, we plot the differences x(T)—\(l0K), 
normalized by the value at 2K. The arrow indicates T c found 
in Ref. [nj. Right panel: initial magnetization curve as a func- 
tion of the applied magnetic field measured at 2K in the 
trigonal phase with a z=0. 546(3) A. The lower critical field 
H c i ~75 Oe is determined from the departure from linearity 
(red line) at low field. The upper critical field H C 2 has been 
estimated at about 7 kOe by the magnetic field at which the 
M-H reverse legs merge at high magnetic field using the crite- 
ria of | AM | < 10 _J emu. The values of H c i and H C 2 include 
the demagnetization effect assumed to be 1/3 in cgs units. 



we infer that a soft phonon mode corresponding to the 
out-of-plane Si vibrations, and very sensitive to buck- 
ling, is responsible for the relatively high value of T c in 
BaSi2 when Si planes get flatter. Moreover, our calcula- 
tions put in evidence a correlation between the softening 
of this phonon mode and the Si-Si bond buckling. We 
expect that such a behavior is a general property of con- 
ventional superconductors presenting corrugated layers. 

Polycrystalline metastable high pressure and high tem- 
perature trigonal samples of BaSi2 were synthesized in 
a Belt-type apparatus from the commercial (Cerac in- 
corporated) orthorhombic phase of BaSi2 (Prima space 
group) with a purity of 98%. Different trigonal struc- 
tures were meta-stabilized by changing the temperature 
and the pressure conditions of synthesis according to 
the phase diagram of BaSi2 d, H3, EH Hfl- The struc- 
tural characterizations were performed by X-ray diffrac- 
tion 9 — 29 measurements done on a Bruker D8 Advance 
powder diffractometer (Kai,2 Cu wavelengths). The Ri- 
etveld analysis (GSAS software) of the X-ray profiles 
confirms the trigonal phase with an average purity cal- 
culated to be more than 98%. The main phase impu- 
rity is the semi-conducting cubic BaSi2 presenting the 
SrSi2-type structure. The three adjustable cell param- 
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FIG. 2: (Color online) Variation of the total energy of trigonal 
BaSi2 per unit cell as a function of the buckling parameter z. 
The low-buckled regime 0.545 < z < 0.575 corresponds to 
values of z that can be explored experimentally. 



eters are the two lattice parameters a and c and the 
z coordinate which represents the 2c?-Wyckoff position 
of Si atoms (z= 0.5 corresponds to completely flat Si 
plans). Among the different samples, the z coordinate 
exhibits sizable variations depending on the pressure and 
on the temperature of synthesis. The two samples pre- 
senting the most pronounced structural difference were 
synthesized at (4.5GPa, 500° C) and (4.5 GPa, 1000° C), 
and have, respectively, [a=4. 061(3) A, c=5. 293(3) A, and 
z=0.565(l)] and [a=4.065(3) A, c=5.347(4)A, and z= 
0.546(3)]. The superconducting transition temperatures 
T c of those samples were measured in a SQUID (QD 
MPMS 5XL). In the left panel of Fig. [fl we can sec 
that the temperature dependence of the zero-field-cooled 
(ZFC) magnetic susceptibility reveals a T c onset of 4K 
for the phase with z=0.565 and 8.7 K for the one with 
z=0.546. For the second sample, with T c close to 9 K the 
M-H loop was measured at 2K (z=0.546), as shown in 
the right panel of Fig. [1] It evidences a type-II supercon- 
ducting state with H c i ~ 75 Oe and H c2 ~ 7kOe. Us- 
ing the Ginzburg-Landau equations flij ] , we estimate the 
penetration depth A ~ 3300 A and the coherence length 
\ = 104 A. 

Our DFT calculations were performed with the ABINIT 
code [20l ] . The exchange-correlation functional was mod- 
eled by a Perdew, Burke, and Erzernhof (PBE) general- 
ized gradient approximation plj . while the electron-ion 
interaction was described by norm-conserving Troullier- 
Martins pseudopotcntials [22j generated with the same 
functional. To obtain the phonon dispersion we em- 
ployed density-functional perturbation theory. Proper 
convergence was ensured using a cut-off energy of 30 Ha 
and a 20 3 /c-mesh with the Monkhorst-Pack sampling of 
the Brillouin zone. A grid of 4 3 q phonon wave vec- 
tors and the tetrahedron technique was used for the in- 
tegration over the Fermi surface. In the trigonal phase 
the lowest energy structure of BaSi2 has lattice param- 
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FIG. 3: (Color online) Panel (a): Fermi surface of BaSi2 in the 
P — 3ml structure [25|]. Panels (b) and (c): two-dimensional 
cuts through the 3D Fermi surface perpendicular to the [001] 
direction for z — 0.55 and 0.58. The increase of the Fermi 
surface when the structure flattens out (smaller z) is evident. 
Panel (d): Electronic density of states as a function of z. 

eters (a = 4.08 A, c = 5.42 A, and z = 0.56 A). These 
values are consistent with our low temperature sample 
within the typical error of PBE calculations. However, 
by inspection of Fig. [2j we realized that the total energy 
curve is rather fiat in an interval of z, analogously to the 
case of CaSi2 [3]. We identify three intervals of values 
of z: a planar region for z < 0.545, a high-buckled for 
z > 0.575, and a low-buckled for 0.545 < z < 0.575. In 
the low-buckled interval, the total energy changes by less 
than 50meV (580 K) per unit cell, explaining why it is 
easy to deform the structure, making it possible to obtain 
experimental samples with different z. 

Despite the flatness of the total energy curve, small 
variations of z in the low-buckled region produce remark- 
able changes in the electronic structure, in the phonon 
modes, and in the EP coupling. The Fermi surface, de- 
picted in Fig. [3Ja) is composed of two parts, a larger 
portion centered at the A point and smaller pockets cen- 
tered at the M-points. The bands composing the former 
have mostly Si 3p xy character, and disperse much more 
strongly in the xy plane than along vertical directions. 
Increasing the buckling of the Si-Si bond leads to a shift 
of these bands, which reduces the size of the Fermi sur- 
face [panels (b) and (c)] and, consequently, decreases the 
value of the density of states at the Fermi level [N(E-p) 
- see panel (d)]. When one reaches z > 0.58, the inter- 
acting orbitals become less confined within the Si layers, 
forming a weak hybridization between the 5d states of 
Ba and the 3p z of Si. Also the pockets, formed by Ba 
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FIG. 4: (Color online) Bottom panel: Eliashberg spectral 
function oi 2 F(uj) (in the inset: EP coupling constant A) as a 
function of z. Top panel: T c in K calculated with the McMil- 
lan formula. The shaded area shows the interval of T c values 
when fi* varies from 0.06 to 0.125. For comparison the dashed 
line represent the experimental values of T c . 



5d and Si 3p xy states, show a light decrease of area upon 
increase of z, but on a smaller scale than the surface cen- 
tered at A. For larger values of z (up to 0.6) the Fermi 
surface does not change significantly from Fig.^c). For 
values of z smaller than 0.55 we can see the appearance 
of a third sheet of the Fermi surface composed purely of 
Si 3s states. This effect can already be seen in Fig. [3jb) 
as a tiny yellow dot centered at the A-point. Note, how- 
ever, that for z < 0.545 the system becomes structurally 
unstable, making therefore this region of limited practi- 
cal importance. Looking at the density of states at the 
Fermi level [Fig. [3jd)] we see that it attains a maximum 
in the low-buckled region (z = 0.545), then it decreases 
with increasing z. These numbers already suggest that, 
from the electronic point of view, superconductivity is 
favored by smaller values of z. 

In order to verify this point, we calculated both the 
phonon modes and the EP coupling constants as a func- 
tion of z. At the equilibrium structure our calculated 
phonon frequencies are close to experimental values [26| . 
For z close to 0.5, i.e. for almost planar Si layers, the 
modes exhibit imaginary frequencies indicating a struc- 
tural instability. In the stable region, most phonon fre- 
quencies show only a weak dependence on buckling. How- 
ever, the A\ g optical phonons, which are mainly com- 
posed of vibrations of the Si atoms along the c axis (buck- 
ling modes), considerably soften upon decreasing z, lead- 
ing eventually to the structural instability. Moreover, the 
higher energy E g mode, composed of Si vibrations in the 
xy planes, is slightly hardened with decreasing z. In spite 
of these changes, the value of rii og , the weighted average 
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of the phonon frequencies, remains basically unchanged 
in the optimal range 0.55 — 0.58, and therefore can not be 
responsible for the variations of T c . At z = 0.56 we obtain 
the maximal value f2i og = 182 K; at z = 0.55 it is slightly 
smaller, f2i og =170K, for z = 0.545 it drops to 156 K and 
for z = 0.57 and z = 0.58 we found f2i og = 179 K and 
171 K, respectively. 

To evaluate the superconducting transition tempera- 
ture T c in the framework of the strong coupling theory of 
superconductivity 27J, \2M we nee d to calculate A, which 
measures the average EP interaction (inset of the bottom 
panel of Fig. @| . As the calculation of A involves averag- 
ing over the Fermi surface, special care has to be taken to 
ensure the convergence of the results. In fact, an insuf- 
ficient sampling can lead to a dramatic underestimation 
of T c 
to 0.1 
rial 



Bai-^Ka-BiOa, etc. This scenario is in full agreement 
with our experimental and theoretical findings of an in- 
crease of T c when Si planes flatten out, and it is compat- 
ible with previous measurements on the disilicide CaSi2 



12j . The screened coulomb interaction fx* was set 



which is a standard choice for this kind of mate- 
The Eliashberg spectral function a 2 F(ui) is shown 
in the bottom panel of Fig. |4j There are two main peaks 
contributing to A: the first, due to the acoustic modes, 
is fairly insensitive to changes in z\ the second peak, due 
to the EP coupling of the A\ g optical modes, moves to 
lower frequency, increasing its area and consequently its 
contribution to A with decreasing z. In the inset of the 
bottom panel of Fig.|4]one can see how this leads to a dra- 
matic increase of the EP coupling constant A with z, and 
therefore to a maximum theoretical T c of around 6 K us- 
ing the McMillian formulation [29[ (top panel of Fig. HJ) - 
Note that the experimental trend for the dependence of 
T c of BaSi2 on the buckling is perfectly reproduced by 
our calculations, which explain also previous experimen- 
tal findings for the similar compound CaSi2- The slight 
underestimation of calculated T c is due to the neglect of 
multi-band effects that are known to enhance supercon- 
ductivity in similar systems, like MgB2 0, 0- 

In conclusion, we presented an experimental and the- 
oretical study of the layered binary silicide BaSi2. We 
show that the buckling of the Si sheets can be modified 
experimentally by using different high-pressure and tem- 
perature conditions of synthesis. The reason for such a 
behavior can be found in our DFT calculations, which 
evidence a broad low-buckling interval where the total 
energy changes by less than 50meV (580 K) per unit cell 
as a function of z. The electronic band structure calcula- 
tions demonstrate that the density of states at the Fermi 
level significantly increases by reducing the buckling of 
the Si planes. However, the flattening of the Si layers is 
limited by a structural instability concomitant with the 
softening of the Si optical buckling phonon mode (Ai s ). 
The coupling of this mode with electrons is also dramati- 
cally enhanced by flattening the Si planes, leading to the 
doubling (from 4 K up to 8.7 K) of the superconducting 
transition temperature. Such competition between su- 
perconductivity and structural distortion can be found in 
a wide variety of conventional superconductors, like the 
Chevrcl-phases, the transition metal carbides, the cubic 



15j,ll6|. Moreover, the mechanism of tuning T c by con- 
trolling the buckling of the layers is likely to be present in 
many other layered superconductors, and therefore can 
provide a new path for optimizing T c . 
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